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УДК 519.642.3 
В.М. ДРАГИЛЕВ 


О РЕДУКЦИИ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ ПЕРВОГО РОДА 
К УРАВНЕНИЮ С ВЫРОЖДЕННЫМ ЯДРОМ 


Рассматривается абстрактная некорректная обратная задача, сведенная К инте- 
гральному уравнению Фредг ольма первого рода. Работа посвящена полуэвристи- 
ческой схеме априорных прогнозов, предназначенной для анализа возможнос тей 
реконструкции по данным об уровне входных погрешности ей. Получены ОЦЕНКИ, 
позволяющие, При подходящем соотношении между параметрами, проводить про- 
гнозы, используя вырожденную аппроксимацию для ядра интегрального операто- 


ра. 
Ключевые слова: некорректные обратные задачи, интегральное уравнение 


Фредгольма первого рода, метод проекции, априорный прогноз. 


Введение. Пусть некорректная обратная задача сведена к интегральному 
уравнению (ИУ) Фредгольма первого рода с вполне непрерывным оператором 
Ь 
[АЧ1© = | (х,5)4(9) 4 =), — хе[с,а], (1) 
а 
где и(х)=и(х)+ди(х), и(х)=[А4 (хх), 4($) е Е а,Ь] - искомая функция 
(оригинал); би(х) - погрешность исходных данных. 


Ввиду некорректности задачи ее практической постановке должен 
предшествовать априорный анализ, состоящий, в идеале, из следующих 
двух этапов. В первую очередь желательно указать некий класс оригина- 
лов, успешная реконструкция в том или ином смысле гарантируется при 
заданном уровне входных погрешностей. Затем, исходя из доступной ин- 
формации об оригинале, следует выяснить его принадлежность к данному 
классу и сделать вывод о возможности (невозможности) успешной рекон- 
струкции. 

Самостоятельный интерес представляет первый этап, который и 
будет в центре внимания. За недостатком априорной информации как пер- 
вый, так и второй этап обычно осуществляется нестрогим образом, с при- 
влечением неких качественных представлений. Простейшим средством 
априорного анализа служат численные эксперименты с имитацией погреш- 
ностей и модельными оригиналами. Другие приемы основаны на исследова- 
нии сингулярных чисел (СЧ) и числа обусловленности [1-4]. По ряду причин 
последний подход может быть наиболее эффективно реализован, когда 
оператор А является конечномерным, т.е. 


баз =б, а =Уф, и, $). г 


т=1 
Подходящая методика априорных прогнозов в основных чертах по- 
казана в статьях [2, 5] и для удобства будет далее изложена. 
Для приложений интересен случай [2], когда ядро ИУ приближает- 
ся к вырожденному, т.е. 
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С(х,5) = С, (х,5) + Ц, (х,5), 
где С,(х,5) - остаточные члены, исчезающие в некотором пределе. 

В процессе априорного анализа заманчиво рассмотреть ИУ с вы- 

рожденным ядром 

Ь 

[А,9](<) = [с, (х,5)4($)4=й<), — хЕ[с,4], (3) 

а 
где, как и в (1), й(х) =и(х)+ди(х), и(х) =[Аа](х), и провести анализ по 
той же схеме, что и в задаче (1), (2). При этом возникает вопрос, насколь- 
ко правомерна редукция ИУ (1) к ИУ (3). С практической точки зрения она 
в какой-то мере уже будет оправдана, если только саму схему прогнозов 
удастся аккуратно распространить на ИУ (3). Цель настоящей работы - 
указать простое условие, при котором такое обобщение становится воз- 
МОЖНЫМ. 

Подчеркнем, что в основе развиваемой схемы лежат качественные 
гипотезы, не имеющие четкой математической формулировки. Естествен- 
но, что и модификацию этой схемы нельзя отнести к чисто математиче- 
ским задачам; в конечном счете речь идет о том, чтобы придать СЧ некую 
разумную интерпретацию. Тем не менее основные соотношения, получен- 
ные ниже математическими средствами, показывают, что условие редукции 
слабо зависит от той неопределенности, которую вносит нестрогий харак- 
тер отправных гипотез. 

Обозначим: А, =А-А,. Тогда имеем: и(х)=и„(х)+ди,(х), где 
ии (Хх) =[Аиа](х), ди,(х) =[А,4](х); оба эти вклада неявно присутствуют 
в правых частях ИУ (1) и (3). Временно пренебрегая вторым вкладом, об- 
ратимся к вспомогательной идеализированной задаче (1), (2), записав ее в 
виде 

[Ам а](х) =им (х)+би(х),  хЕ[с, 4]; (4) 
известной здесь считается только сумма й(х)=и,(х)+ди(х). 
Сводка прежних результатов, относящихся к задаче (4). Выбрав ор- 
тонормированный базис {},(5)}” в пространстве Г,[а,6], обобщенные 


решения ИУ (4) можно строить проекционным методом в виде 
4, =У`Е1,(), (М<М). (5) 


Здесь коэффициенты &, отыскиваются методом наименьших квадратов, 
т.е. минимизируют невязку ИУ в пространстве [,[с,4]. В последующем 


контексте указанная процедура реконструкции может считаться чисто во- 
ображаемой (априорные прогнозы сводятся к расчету СЧ и их интерпрета- 


ции). Подразумевается, что системы {ф„(х)}"_ и {и „(5)}"_, обладают ли- 
нейной независимостью, а проекции функций [,(5), (п=1,2,...,М) в ли- 
нейную оболочку функций и,„(5), (т=1,2,...,М) образуют базис этого 
подпространства. 

Обозначим: Н„, - подпространство, натянутое на систему 


И Ну - его ортогональное дополнение до Г.[а,6]; Р и 


м 


294 


Вестник ДГТУ, 2006. Т.6. №4(31) 155М 1992-5980 








Р,=1-Р, - соответственно проекторы из Г.[а,] в Н, и Н,, 
Ви =(„А, АЛ.) ‚ Вх - эрмитова матрица с элементами (1<т,п< М); 
с” - СЧ матрицы Ву, упорядоченные по убыванию: с" > о’ 2..2 оу 
(можно показать, что в наших условиях с” >0); С, =д"“/оу - число 
обусловленности матрицы Ву„; 4,(5)=Р,а(5); а„(5) - решение (5) за- 
дачи (4) с би=0 в правой части. Отметим, что 4’(5)=а,(5)+Р.Р, 4(5), 
где О, - некий линейный оператор [5]. Для функций, определяемых на 
отрезках [а,] или [с,4], символы ||| будут означать соответственно 
норму в Г, [а,6] или Ё.[с,4]. 
Введем в рассмотрение относительные погрешности: 
д= |6, [| п» =|а-а,|Иа[, п» = 9 - а Иа о = | - 4/91, 
Я» = |4» —ч|/ 
Погрешность решения (5) й, можно оценить при априорно задан- 
ных погрешностях 77, и Д. В статье [5] получены оценки: 
Йы < Пу + Пу +} (6) 
У (7) 
где д. =|р,|. 
Для погрешности 77" имеется оценка [6, 7] 

у (8) 
где у/=|м, |/|м.|, и, =А,а,. Фактор у, трактуется как функционал от 
4(5). Отметим, что у, =ТГ и у >1 при №М<М [7]. Число обусловленности 
С, монотонно растет вместе с № в силу принципа промежуточности [8]. 

В канонической версии метода проекций базис образует собствен- 
ные функции (СФ) эрмитова оператора А» А„, что делает метод регуля- 
ризующим [9] и дает оптимальные значения у, =0, #7" =п, [5]. При ином 
выборе базиса, вообще говоря, 77, >17], , т.е. возникает особого рода доба- 
вочная погрешность. В связи с этим предполагается, что для всех М < М 
значения нормы у, меньше или порядка единицы. Согласно (7) при таком 
условии погрешности 77, и 7», всегда сопоставимы, и тем самым эффек- 


тивность проекционного метода сравнима с эффективностью регуляризу- 
ющих схем. Волевой выбор базиса позволяет обойтись без исследования 
СФ, что особенно привлекательно в многопараметрических задачах. По- 
путно поясним смысл редукции ИУ (1) к ИУ (3): в канонической версии ме- 
тода проекций она упрощает отыскание СЧ и СФ, а в обобщенной версии - 
оценивание нормы х, [5]. 


При помощи оценок (6)-(8) удается анализировать возможности 
реконструкции в рамках следующих упрощающих гипотез. 
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1) Предлагается пренебречь множителем у, в оценке (8) [2, 4]. 
Эвристическим основанием для этого служит тот факт, что функционал 
у, обращается в единицу на всех ненулевых оригиналах 4(5)еЕН, и 
непрерывен на каждом из таких элементов; вместе с тем при любой раз- 
мерности решения № успешно восстанавливаться могут только оригина- 
лы, близкие к функциям класса 4(5)=Н,. 


2) При С, >>1 и при случайном характере погрешностей ди(х) 


оценка (8) чаще всего будет многократно завышена даже без учета множи- 
теля у,. Утверждается, что максимальная размерность №, при которой 


и 


погрешность 7” обычно еще остается много меньшей единицы, может 
быть найдена по критерию [1, 2, 4] 
М = тах {№: СА. (9) 


тах 
Высказывание о том, что погрешность 7” скорее всего мала, явля- 


ется чисто качественным; получить оценку более точную, чем (8), можно 
только при наличии дополнительной априорной информации. 

Согласно гипотезам 1) и 2) успешную реконструкцию допускают, 
как правило, те и только те оригиналы 4(5), которые при М№М=М№„„ адек- 
ватно приближаются своей проекцией 4,(5); последнее предполагает, что 


77, <<1. Учитывая оценки (6), (7), при М=мМ 


х 


для этих оригиналов мож- 


тах 
но прогнозировать вероятное восстановление с погрешностью й, много 
меньшей единицы. При такой погрешности качество изображения 4, (5) 
будет примерно соответствовать качеству аппроксимации 4,(5). 


Тем самым создается возможность делать прогнозы, исходя из за- 
данного уровня входных погрешностей Д и вычисляемого конечного набо- 
ра величин - значений числа обусловленности. 

Оценка погрешности решения в задаче (3). Пусть теперь (5) есть по- 
строенное тем же способом обобщенное решение ИУ (3). Последнее отли- 
чается от ИУ (4) только дополнительным вкладом би,(х) в правой части. 


Рассматривая этот вклад как добавку к функции би(х) и применив нера- 
венство треугольника |6и + ди, | < [би| + |би, 








‚ при помощи оценок (6)-(8) 


получим такую же систему оценок с тем отличием, что в правой части не- 
равенства (6) добавляется слагаемое й, =(+7,)7,С,А,, 


где А, =|би, 








Им,|. Функция 4,(5), как и прежде, минимизирует невязку 
ИУ (4) при ди=0. Из ее сингулярного разложения найдем, что 




















оч. |<|м„|. Кроме того, имеем |0и,|<|А/|4|. Следовательно, 
А. <С,(|АИ а )-(|4|/.|). Здесь, по принципу промежуточности, 
о" >о0'. Обозначим: Д,=|а|/|а.|, а’>|А,| - произвольная оценка 


сверху для нормы |А, 








‚ "=, а=а’/ ам’. (Отметим, что число 
а" =, А,А,Х), есть оценка снизу для нормы |А,|). 
Из предыдущего вытекает набор неравенств 
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д <арВ.С., (№М=12....М ). (10) 


Введем дополнительное обозначение АД = /а . Полагая в (10) 
№=М№„„„ И учитывая, что при этом АС, < 1, получим: 





А, <А* В, С, <А*В,/А= В. (АГАРА, (М=М). (11) 
Наконец, используя (11), приходим к оценке 
п, <(1+1.)у, В, (А!АУ, (М=Мщ). (12) 


Интерпретация оценок (11), (12). В рамках принятых допущений пер- 
вые три сомножителя в (12) сопоставимы с единицей для всех оригиналов, 
для которых может прогнозироваться успешная реконструкция. Тем самым 


при соотношении Д >> А вся схема прогнозов дословно переносится на ИУ 
(3) с тем отличием, что в оценке для Й, (6) появляется добавочный член 
7, <<1. Практически достаточно, чтобы отношение А/А было больше 
10 ; тогда предполагаемые значения 7, будут составлять порядка 0.01 и 


менее, что вполне отвечает качественным представлениям о малой по- 
грешности в духе гипотезы 2). 

Сделанный вывод нуждается, правда, в одном уточнении. Дело в 
том, что относительная погрешность исходных данных до сих пор была 


определена как л=|би |/|м, |, в реальности же приходится иметь дело с 
априорно заданной величиной л=|би|/|и|. Это обстоятельство будет 
несущественным, если только, согласно априорной информации, 
[би | << |и|. В самом деле, при соотношениях Д>>Д, 2, <<1 из (11) 
находим, что |ди, | <<|6и|; отсюда при |5и|<<|и| следует, что 


|= [и |. Небольшие различия в определении А могут влиять на вычис- 


[мы 
ляемую по критерию (9) размерность №. Это, однако, не играет роли, 


так как сам критерий является приблизительным, а критическое значение 
С,А=1 - отчасти условным. Применение критерия становится наиболее 


содержательным, когда число обусловленности быстро возрастает, т.е. 
С. >>С, В некотором интервале значений №, в пределах которого 
находится №. [4]. В таких условиях даже большие (скажем, на +50 %) 
вариации А способны изменить значение № „, не более чем на +1. 
Неопределенность, возникающую при формулировке критерия (9), 
можно выразить математически, заменив условие С,4<1 на С,/А<п, где 


"- константа, регулирующая максимальные значения погрешности 7”. 
Тогда неравенства (11), (12) останутся верными с точностью до замены 
д=\Уа на Д= 7“ и добавления множителя 77 в правой части (12). 
При 7<1 условие АД >> ДА только лишь ослабляется. Задавать же большие 
значения 7>>1 на практике не имело бы смысла: при таком выборе апри- 


орные оценки не могут гарантировать даже отдаленного качественного 
сходства между решением (5) и проекцией 4,(5). При отказе от гипотезы 


1), т.е. при учете фактора у„,, возникнут конкретные числовые ограниче- 
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ния на погрешность 7„,, при которых значения у, еще остаются близкими 
к единице. Это, однако, никак не отражается на применимости условия 


А>>А. 
Анализ частного случая. Ясно, что редукция ИУ (1) к ИУ (3) в целом 
ухудшает возможности реконструкции; результаты настоящей работы не 
позволяют строгим образом контролировать эти последствия. Укажем все 
же способ для косвенного контроля, предварительно рассмотрев следую- 
щий пример. 

Пусть {4 (5)}” есть ортонормированный набор СФ оператора 


А”А, отвечающих его СЧ д, > о, >..., Р,, - проектор в линейную оболочку 
СФ {Р (5). Предположим, что би, /би << 9/01, т.е. в спектре СЧ 
имеется резкий скачок при т = М . Выберем аппроксимацию интегрально- 
го оператора в виде А,=АР,; для этого в формуле (2) положим 
и,(5)=Ф (5), 6, (х)=[АЧ, (>). В качестве базиса {},(5)}"_, выберем 


систему СФ {Ч,(5)}"_ . В этих условиях имеем: |А„|=о., |А, 








= Ома, 


Си=о,/ ви, @-|А,|. Заметим, что /а>\/о„„/о,, или, что равно- 


сильно, Д>Д’/Е, где А’=о,,/ои, ==00,./ 0 <<1. Условие 


А>> А выделяет область погрешностей А>> Д'. Легко видеть, что в этой 
области при любом методе решения ИУ (1) исключается успешное восста- 
новление каких-либо иных оригиналов, кроме функций двух типов: таких, у 


которых в разложении по базису {Ч (5)}” преобладают члены с т<М, 


и таких, у которых все эти члены подавлены. Функции первого типа допус- 
кают реконструкцию на основе ИУ (3) и улавливаются при априорных про- 
гнозах. При решении ИУ (1) методом А.Н. Тихонова [10] возможна рекон- 
струкция также и функций второго типа. Но если, при всех прежних усло- 
виях, система СФ подходящим образом упорядочена (например, по возрас- 
тающему числу нулей или экстремумов), то функции второго типа не будут 
представлять большого прикладного интереса. 


Далее для простоты положим а’ =|А, 








: Тогда 
Ао ЕС, так что Д<<1/С,. При этом условия 


А<<А<1/С, совместны, и в таком диапазоне погрешностей критерий 
(9) дает М, =М. 
В общем же случае, т.е. когда А, = АР’ , условие А>>А приво- 


т: 


дит к ограничению №„<М, где М=шах{М№: С, А<П. Если окажется, 


1<№М<М 
что №М<М, то это может означать, что аппроксимация интегрального 
оператора неудачна (либо выбранная оценка с’для нормы |А, недоста- 








точно точна). Если ж  М=М, то тогда можно надеяться, что при 
А>>А вся ситуация сходна с рассмотренным примером. А именно: редук- 
ция ИУ исключает восстановление функций класса 4(5)=Н, и близких по 
норме, но не оказывает большого влияния на восстановление функций 
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класса 4(5)=Н,„, и близких по норме; реконструкция же всех остальных 


функций объективно затруднена. 

Замечания. Фигурировавший в рассуждениях проекционный алгоритм с 
выбором № по критерию (9) не столь эффективен, как метод А.Н. Тихоно- 
ва с апостериорным выбором параметра регуляризации (тем более, что 
последний метод не требует и редукции ИУ). Использовать такой алгоритм 
непосредственно для решения обратной задачи целесообразно только 
лишь в специальных случаях при наличии подходящей априорной инфор- 
мации. Тем не менее, априорные прогнозы дают представление также и о 
возможностях других методов применительно к восстановлению функций 
класса 49(5)=Н, и близких по норме. Это высказывание подкрепляется как 


модельными примерами [21, так и качественными соображениями. Его 
краткая аргументация сводится к тому, что даже при волевом выборе ба- 


зиса СЧ о” имеют обычный смысл: они выступают весовыми коэффициен- 
тами, с которыми вклады от ортогональных проекций (} 


„Чл, 1,(5) 
(п=1,2,..„М№) входят в правую часть ИУ (4). 

Если исходные данные и(х) задаются в виде / - мерного массива 
значений в неких опорных точках х,, то тогда операторы А, А„ удобно 


рассматривать как действующие из [..[а,в] в Ю” , соответствующим обра- 


зом переопределив нормы и операцию эрмитова сопряжения. 
Выводы. Условие, делающее целесообразным редукцию интегрального 
уравнения первого рода к уравнению с вырожденным ядром, получено в 
форме ограничения снизу на уровень входных погрешностей. При таком 
условии класс оригиналов, успешное восстановление которых может про- 
гнозироваться в рамках рассмотренной схемы на основе редукции, близок 
к аналогичному классу оригиналов, допускающих успешную реконструкцию 
в идеализированной обратной задаче с вырожденным ядром. 
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\У.М. ОВАСТЕЕУ 


ОМ ТНЕ 1ТМТЕСКАЕ ЕОЧАТТОМ ОЕТНЕ ЕТК$Т КТМО 
ВЕООСТТОМ ТО ТНЕ ЕОЧАТТОМ МТТН А ОЕСЕМЕКАТЕ КЕКМЕЕ 


Ап абзгасЕ -розе туегзе ргоБет гедисе ю ше тедга! Егеддпойт едчца- 
Чоп оЁ Че Яг$Е Кта 15 сопзаегед. ТВе рарег 15 деус{еа тю {Пе зет!-Веий$Ис 
эспете оф {Не а рйой рге4своп$ аезЧптеа Гог {Пе апа!уз!5 оЁ те гесоп&гисНоп 
ро Без гот {е д№еп |е\е| о {Пе па! даа еггог. ТПе иррег Боипа$ аге 
деймуе у/ИсИ аатй, п {Ие сазе оЁ Ме арргорйае геаНоп Бебмееп {пе ра- 
гатеегз, {© саггу оц {Ле рге-ЧсНоп$ изтд {Ле дедепегае арргохйтаНоп Юг 
{пе Кегпе! ог 1е Итедга! орегаюг. 
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